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Abstract 

We use recent results regarding the geometry of the U(3) bosonic algebraic model to 
motivate a hybrid algebraic-Schrodinger approach to the 2D Franck-Condon problem anal- 
ogous to ID approaches. This approach allows one to analyze bent to linear Franck-Condon 
transitions of triatomic molecules as well as clarifies the angular momentum dependance of 
the Franck-Condon intensities. 



1 Introduction 

The study of transition amplitudes for molecules with initial and final states given by differ- 
ent electronic configurations has held great interest. In the Franck-Condon limit one models 
the two configurations by similar effective potentials with different geometric parameters. 
The transition intensities are then given by the wavefunction overlap of the two potentials 
(sudden approximation) |J. 

Recently, the Franck-Condon (FC) problem of polyatomic molecules has been made 
tractable by the introduction of a hybrid algebraic-Schrodinger approach In this 

approach a spectrum generating algebra of Ui(2) x ■ • • x £4(2) (where there are k bonds) is 
used to obtain wavefunctions || 

\4>)= E c^Jfivj,^)®-.-®!^],^), (i.i) 
»i, •••>** 

by fitting spectra of a molecule in a particular configuration. In general for the two different 
configurations the zth bond contribution could change from the representation [iVj] to [N-]. 
The Franck-Condon factors for the polyatomic case are realized in terms of the single bond 
factors by calculating the matrix element {ip'\U\%lj) |2|||, where 

tr = f(^,A 1 )®...®f(^,A fc ) (1.2) 

and the operator t is defined by its matrix elements as determined by Schrodinger overlaps: 

([N'],n'\t(*,A)\[N},n) = / dx^'(a'; x)^(a;x - A). (1.3) 
a' J 
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The parameters ^ and A are typically fit from transition spectra. Due to the U(2) D 0(2) 
chain's correspondence with a Morse oscillator || one often uses Morse eigenfunctions with 
anharmonicity parameter ^-j- for Schrodinger wavef unctions. It is worth noting that the 
heuristic approximation in || showed that this can be equivalently thought of as simple 
harmonic oscillator (SHO) overlaps where the scale parameters acquire a linear n dependent 
correction with coefficient 0(1/N). That is, one would use the equation for the overlap of 
SHO wavefunctions, which has no N dependence, and substitute a scale a = a (l— (,(n+^)). 
«o was taken as the scale of the SHO best approximating the Morse and £ is 0(1/N). In 
this way one introduces geometric parameters not contained in the algebra. 

Although this approach has been successful it necessitates that the initial and final 
configurations not be too dissimilar — in particular they must have the same normal mode 
expansion. This precludes, for example, a transition between a linear and bent triatomic 
molecule (Figure 1). In the linear configuration the center atom has two normal modes. In 
the bent configuration one of the normal mode solutions is spurious — corresponding to an 
overall rotational degree of freedom. Because of this, the treatment of vibrational modes 
in each configuration necessitates the freezing of a different number of degrees of freedom, 
whereas when studying transitions between each configuration one needs to have the same 
degrees active. 

One may attack this problem by using more complicated algebraic models which have 
both geometries built in. In this way some parameters which before were artificially inserted 
through the hybrid method become natural — i.e. since they are built in the algebra they 
are determined by spectra and do not need to be fit with transition data. 

However, in these more complicated algebraic models, exact Schrodinger correspon- 
dences don't in general exist. Thus, this approach necessitates an alternate interpretation 
of algebraic parameters as geometric, configuration space quantities. Recently such an in- 
terpretation has been provided |7[] by exploiting the many approximate correspondences 
between algebraic and Schrodinger pictures. We use these results to develop a hybrid ap- 
proach for the bent to linear transitions in Figure 1. This is most easily done by using a 
spectrum generating algebra of Ui(2) x Z7 (3) x U 2 (2) (a U(2) for each bond and the U(3) 
for the two additional degrees of freedom of the center atom). The FC transition operator 



would then be U — t\ ® T ® i 2 , where ti, t 2 are defined by equation [L3| . It is the purpose 
of this publication to motivate the definition of T (eqn. p.5| ). 



We begin by demonstrating, via a coherent state limit analysis 0, that U(3) is an ap- 
propriate algebra to describe the degrees of freedom of the center atom of the triatomic. 
We proceed with numerical studies to study the implications and discrepancies of this in- 
terpertation. Finally, we reconcile these discrepancies using [F| which naturally motivates 
our definition of T. 



2 The Franck- Condon Problem and U(3) 
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2.1 Overview of U(S) Statements on 2D Problem 



Th algebraic approach for 2D problems was presented by [Q. One considers symmetric 
(bosonic) representations of u(3). There are two chains of interest: 

U(3)dU(2)dO(2) I 
U(3) D 0(3) D 0(2) II ' 

We use the same notation as |§ for the generators except chose a different 0(3) subgroup 
as explained in J7|. Please note that the 0(3) group is a dynamical symmetry subgroup and 
does not have the interpretation of a rotation in configuration space. 
The general Hamiltonian of the U(3) model is 

H = en + 5n(h + l)+(3l 2 -AW 2 , (2.2) 



where e, 5, and A are taken as positive or 0. Setting A = (e = 5 = 0) gives a Hamiltonian 
with dynamical symmetry I (II). The spectra for each dynamic symmetry may be determined 
from the well known solution to the branching problem for a symmetric representation of 
u(3) labelled by N (the eigenvalue of the u(3) Casimir n + n s ) ||. The basis corresponding 
to chain I is labeled by the eigenvalues of the u(2) and o(2) Casimirs n and I respectively. 
The basis corresponding to chain II is labeled by the o(3) Casimir's eigenvalues, uj(oj + 1), 
and again by I. 

The spectra of each chain led the authors of || to the interpretation of each dynamical 
symmetry as an azimuthally symmetric potential with minimum at radius (chain I) and 
at non-zero radius (chain II). 



2.2 U(3) Coherent State Limit 

The interpretation of || is reaffirmed by simply studying the classical coherent state limit 
calculated in 0. 

Taking the coherent state limit of the Hamiltonian [2.2| , setting all momenta to zero and 
dropping additive constants one finds the potential in group coordinates JTOf to be (up to a 
meaningless multiplicative factor): 

V c \(r) = r, l -r 2 + l -r A (2.3) 

where 

e + 25 + P-4A(N-l) 
1 (5 + AA)(N-l) ' y ' } 

It is easy to compute the position of the potential minima 

°r~ "'I (2-5) 

1—1] 7] < 
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Thus we find: 



^ci(r min ) = { ^ ° • (2.6) 



That is we have a second-order phase transition at rj = 0, or equivalently 4A(N — 1] 



e + 25 + p [11 



The exact same analysis may be carried out in projective coordinates [TO] revealing: 



. .. 7] > 1 

Vd'( r min) = i 1 (l-,,) 2 < 1 , (2.7) 




1+J?' 



where now 

= (e + 2£ + g) , = _5_ 
77 ~ 4A(iV - 1) ' 77 ~ 4A 

and 



(2. 



e,„= a (2.9) 

Thus we again find a second-order phase transition at AA(N — 1) = e + 25 + /3. 

In either coordinates the potential minima moves from r = (corresponding to a linear 
configuration) to r ^ (corresponding to a bent configuration) at the critical point. We 
conclude that the algebraic model is rich enough to include both geometries depicted in 
Figure 1. Algebraic Hamiltonians having nearly a U{2) dynamical symmetry correspond to 
the center atom in a linear triatomic, whereas those near the 0(3) limit correspond to the 
center atom in a bent triatomic 



2.3 U(3) and Schrodinger FC Connections 

The FC factor for the two configurations is easily studied from a Schrodinger perspective. 
Assuming that, whatever the actual nature of the potentials, they may be approximated 
about their minima as harmonic we may numerically calculate the FC overlaps. The results 
of such a calculation are displayed in Figure 2. The parameters (Table 2) are chosen to be 
relevant to the bent to linear FC transition ( 1 I?2 — > £+) of CS2. The frequency of the linear 
configuration was taken from reference ||12|| . The harmonic distance scale, was deduced 
by assuming the effective mass was that of the carbon nucleus. The radial displacement 
was deduced from the geometry of the bent configuration as published in reference |T3"| 
assuming that the heavy (S — S) axis was essentially stationary. The associated distance 
scale of the bent configuration was assumed to be the same as the linear scale. Since CS 2 is a 
particularly shallow molecule we have included an additional plot (Figure 3) to demonstrate 
the behavior for a larger radial displacement. 

If one models the potential as exactly harmonic the displaced oscillator potential has 
the idiosyncrasy of a 'cusp' at r = 0. This is of little concern since the effective poten- 
tial is dominated by the angular momentum barrier at this point. Approximate analytic 
expressions may be obtained for this limit as detailed in Appendix A. 
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We wish to emphasize that although these graphs serve as a valid starting point for a 
more complete analysis of CS2 transition intensities their primary purpose here is heuristic. 
A full analysis requires a careful fitting of the bent configuration distance scale. Additionally, 
the bending modes considered here are known to couple strongly to symmetric stretching 



modes |T2|] — i.e. a full analysis would require coupling additional C/(2)'s as discussed in the 
introduction. 

One may consider more realistic bending potentials such as a Poschl- Teller: 

V = V [l - cosh" 2 a(r - r*)] . (2.10) 

In this case the cusp at the origin still exists for the displaced oscillator (non-zero r*) but 
is tamed due to the long distance flattening of the potential. In the limit where the minima 
is far from the origin the cusp essentially vanishes. Figure 3 shows numerical results for a 
Poschl- Teller model of a bent to linear transition. The parameters are chosen such that the 
the potentials are approximated to second order by exactly the harmonic plots included in 
the same figure. That is, the harmonic distance scale, a 4 = 2mVod 2 /h 2 is set to the same 

2 

value as the SHO FC factors. The remaining parameter, taken unitless measure 

of well depth. It was chosen to be sufficiently small to emphasize differences between the 
FC factors of the two potentials. 

The previous section's analysis implies that the FC factors for a bent to linear con- 
figuration in the algebraic picture are given by exactly the inner product of the algebraic 
wavefunctions for hamiltonians near the 0(3) chain (bent configuration) and on the U(2) 
chain (linear configuration). The overlaps for several such 'bent' hamiltonians are given in 
Figure 4. The algebraic 'bent' hamiltonian was taken to be of the form 

H=(l-0n-j^- T) W 2 . (2.11) 

The parameter £ was chosen to match the intensity maximum with that of several harmonic 
Schrodinger calculations (corresponding to unitless radial displacements of 3, 5, and 7). The 
results of calculations for two significantly different irreps., [N], are shown to emphasize that 
the structure is generic and not a function of any special choice of parameters. 

Comparing Figures 2-4 again reaffirms the interpretation of the two U (3) chains as bent 
and linear configurations of a 2D problem. Comparing the SHO and Poschl- Teller figures 
one notes although not in exact agreement they are very similar given the large differences 
of the Schrodinger potentials. Qualitatively the £7(3) graphs also appear similar with the 
possible exceptions of (1) their (expected) truncation at higher n; (2) their dramatically 
sharper peaks than the Schrodinger FC graphs; (3) their amplitude's diminished sensitivity 
to the amount of radial displacement. 

2.4 Scale Changes 

The numerical study raises two questions (1) What are the relations between the algebraic 
parameters (determining the hamiltonian's proximity to either chain) and the Schrodingers 
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(more geometric) parameters? (2) Is there an interpretation for the apparent qualitative 
differences (the sharper peak) between the algebraic and Schrodinger pictures? 

Both these questions can be addressed by considering the results of 0. The two relevant 
results, which we reproduce here, include the intrinsic distance scale of the harmonic approx- 
imation to a Hamiltonian near the 0(3) limit and an expression for the radial displacement 
for the same Hamiltonian: 



a 2 ^2 m "° (3) 




where X = e + 25 + (3, Y 



m(jJ (3) 
5(N - 1), and Z 



(2.12) 



(2.13) 



A A (N — 1), and the condition that we 



are near the 0(3) limit implies ^ is small. The parameter muj has the interpretation of the 
ratio of distance to momenta scales, i.e. WZJ = a(mV )2 for a potential V = Vof(ax). 

Additionally the results of |7j imply that Hamiltonians from either dynamical symmetry 
not only correspond to different geometries as implied by |2.2| but additionaly to different 
intrinsic scales (see Appendix |C]): 



Cn 



rrvjJu{2) 2 log 2 



mu (3) 2 - log 2 



4 f31ogiV . , , . r 

1 + log2(2-log2) Uir + (nl0gl ° g2 + C) N 



c = log log 2 — log 



2f 



2-log2 



(2.14) 
(2.15) 



This implies that the algebraic overlap of a U(2) Hamiltonian with an 0(3) one is not simply 
analogous to the overlap of radially displaced oscillators, but analagous to the matrix ele- 
ments of an operator which radially displaces and dilatates (much like the operator matrix 
elements calculated in ||) changing the natural scale of the problem. The degree of dilata- 
tion depends on the proximity of the second hamiltonian to either chain. For chains near 
£/(2) the dilatation paramater is essentially 1. As one moves nearer to the 0(3) chain the 
dilatation parameter increases, approaching the value given by equation [2.14| . This effect 
must be accounted for in any algebraic or hybrid approach to the FC problem. 

We saw in the introduction that in the ID problem the scale parameters may have been 
thought as harmonic scale parameters with corrections linear in the quantum number n of 
order 1/N. The scenario is similar here — except the quantum number is now n and there 
are additional corrections of the larger order log N/N. 

2.5 Schrodinger and Algebraic Parameter Relations 



Equations [2.12j and [2.13| for the harmonic dilatation and radial displacement establish the 
needed connection between algebraic and geometrical parameters — at least in the regime 
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where we are near the 0(3) limit. These quantities can be easily related to experimental 
data. 

Experimentally one can find the lower energy level spacing AE cxp , the reduced mass of 
the particle m exp , and from rotational spectra the distance of displacement r exp . In terms 
of these quantities one may compute the unitless distance using the harmonic oscillator 
dilatation 



< = 6XP .l 6XP - (2-16) 



m„ p AK 
It 

Equating al xp r^ xp with the previous expressions (a 2 r* 2 given by |2.12| and |2.13| ) one finds 



2N [1 + corrections] = ^^ r^. (2.17) 



This could be very valuable when fitting spectra. Since the 0(3) chain represents the 
'maximum' radial displacement @ this expression gives a lower bound for N. One may 
begin fitting data for the N which satisfies this equation with the corrections set equal to 
0. If the 0(3) chain spectra with this N doesn't fit the experimental data then one can try 
higher N and move off the 0(3) chain a corresponding amount so that the equation with 
the corrections is still satisfied. 

3 2-D Franck- Condon Problem: The Prescription 
3.1 Adding Dilatations 

We are now ready to develop a prescription for calculating the FC factors for a U(3) algebraic 
model. We begin by considering the FC problem of two configurations of a molecule both 
described by the Z7" (2) chain — i.e. two linear triatomics. Following the ID procedure we 
propose a hybrid approach based upon calculating the matrix elements of the operator T 
defined in terms of 2D SHO Schrodinger overlaps: 

([N'),n',l\f\[N],n,l) =T n , <l (- / ) (3.1) 

a 



The appropriate overlaps are calculated in Appendix [D] (equations |U.3| and |D.2| ). T by 
construction does not connect subspaces of different I. 

3.2 Final Procedure 

We have yet to add any dependance on the representation label N. One could add such 
dependance by hand, appealing to an analogy with the ID case. However, when one has a 



bent to linear transition such an appeal is largely unnecessary due to equation 2.14 . 

For simplicity suppose we have a molecule with a bent configuration whose spectra is fit 
by a hamiltonian of the 0(3) chain in representation N and a linear configuration whose 
spectra is fit by the U(2) dynamical symmetry with 5 = (3 = (in this scenario the second 
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representation label N' can be arbitrary since the hamiltonian will generate identical low 
lying (experimentally measurable) spectra regardless of N'). Further, suppose the induced 
harmonic dilatations a and a' (given by equation 2.16| ) can be calculated. In this instance 



we know that expanding the 0(3) basis in terms of the U{2) basis is equivalent to expanding 
in terms of a SHO with harmonic dilatation a^ 2 ) = |Cn« 2 (where ( n is determined by [2.14Q . 
Thus the FC transitions should be the matrix elements: 

([N'],n',l\f\[N},u;,l), (3.2) 

where T is defined by 



([N% n', l\f\ [N], n, I) = T Ui<l \-J -( n . (3.3) 




Expanding in large N we see J\Cn = a + bn from equation 2.14 where b is O(-k) and a 



2 

has constant contributions and ^ff^- contributions. With the exception that a ^ 1 this is 
exactly the correction one has in the ID case. That is, we have just let a — > a(a + bn). 

In this extreme case the FC factors can be calculated with no extra fitting parameters. 
However, our calculations depended on expansions in large N and thus we would expect 
them to be inaccurate for larger n. This can be compensated for phenomelogically by 
allowing b to be fit to compensate for ignored terms. 

Although the above situation is only for a limiting case, these results (along with insight 
from the ID analysis) will imply exactly what will happen in other more realistic cases. 

• Bent configuration near the 0(3) basis: One can in principle calculate the scale 
dependance near 0(3): rTKU|ofro(3) ~ ™^|o(3)(l + + /i~) where v and \i would 
have to be determined by repeating the calculation of |7| to second order. In such 
a calculation, it is clear that the result will go like a — > a(a + bn + dl 2 ), where 
a = a(l + 7) and 7 is 0(^, b — 6(1 + 7') and 7' is 0(^, -%); the new parameter d, 
introducing I dependence, is 0(-^a, 7^2); and we have ignored a n-l cross term which 
is of a significantly smaller order. Although these coefficients could be calculated in 
principle, one must phenomenologically fit b and d in order to describe higher states 
anyways. Since the I dependance is of lower order it seems likely that d may not even 
be calculable with todays data. 

• Linear configuration in U(2) basis but 5^0: In this instance the u(2) requan- 
tization was approximate and inadequate for higher eigenstates. However, compar- 
ing the spectra to a Dunham expansion one sees that - plays the role of an anhar- 
monicity parameter (recall equation |2.2|). Thus, from the ID case one concludes that 



a' — > a'(l + b'(n + 1)) where b' must be fit but should be of the order of 



s 



Either configuration far from both chains: Although our scheme has not allowed 
us to do explicit calculations in this regime, in terms of phenomologically fit parameters 
the result should still be clear. One should replace the induced harmonic dilatation 



S 



by a — > a(a + b(n + 1) + dl 2 ). In this scenario a should be between 1 and a depending 
on the proximity to either chain, b has contributions due to both 5 anharmonicities 
and the 0(3) chain, and one again expects that d is of significantly smaller order. In 
this regime all parameters must be fit to data. 

Let us recapitulate. The scheme we propose involves fitting the spectra of each configu- 
ration of the molecule in a U(3) model to obtain wavefunctions \[N],i/)e), \[N'],iP'e') ■ The 
FC factors are then 

([N'W EI \f\[NU E ) (3.4) 
where we define T in terms of the U(2) basis: 

an mm*, o = * v t„ „ , ( ) . (3. 5) 

U J ' ' ' \a'(d' + b'(n' + l) + d'l 2 )J V ; 

In general the parameters must be fit, but whenever any of the limiting cases appear the 
appropriate theoretical values may be substituted. For example, if the linear configuration is 
highly harmonic one has a' = 1 and b' = d' = 0. Further, since we expect the I 2 dependance 
to be nearly negligible and we have a I = I' selection rule, we may expand and combine d 
and d' into one parameter. 

Note that this scheme may also be used for bent to bent transitions. This would be 
more useful than the ID procedure @,[|] if one was interested in the I dependance of the 
transitions for instance. 



4 Conclusions 

Using the result of [0] involving the geometry and scales of the dynamical symmetries of 
U(3) we have explained the difference one sees when considering FC overlaps as described by 
the U (3) algebra and FC overlaps as described by radially displaced Schrodinger oscillators. 
Additionally, the analysis has given us the practical result of a minimum value of N as a 
function of simple experimentally measured quantities. 

We observed that this scale dependency leads very naturally to a description of two di- 
mensional FC factors. To compute the FC factors one expands the algebraic wavefunction in 
the U(2) basis and formally replaces the basis element's overlap by the 2D SHO Schrodinger 
overlap. Corrections for N, anharmonicities, and even I dependance are made by letting 
the SHO dilatation constant have n and I 2 dependent contributions. The resulting formulas 
are analogous to the ID results obtained by expanding in anharmonicities. 
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Appendices 

A Approximate Analytic Expressions for 2D FC over- 
laps 

A.l Linear Configuration 

In the linear configuration in the harmonic limit one has, after separating out the azimuthal 
portion of the wavefunction, the equation: 

-4- ( r4 v ( r ) \ ~ +(E- X 2 r 2 ) v(r) = 0, (A.l) 

r dr \ dr J r l v ' 

where A is related to the frequency of the oscillator by A = ^ and E = j£E. Subject to 
the boundary conditions v(r)\ = 0, rv(r)\ r=Q = the solution is well known to be fbj : 



with eigenvalue 



\ (\l\+n r )\ " 6 L ^ Xr h (A.2) 



E = 2A(|/| + 1 + 2n r ). (A.3) 



A. 2 Bent Configuration 

For a system w/ equilibrium located at some r = r we again assume that for the lowest 
states the system is well approximated by a harmonic oscillator. Our wave equation is thus: 

\Tr [ r Tr V{r) ) ~ ^ V{r) + ( & ~ X ' 2{r ~ ro)2 ) V{r) = °- (A ' 4) 

Notice, that this mild transformation will greatly change the form of the solutions (as 
opposed to the ID case) since the laplacian is not invariant under radial displacements. Note 
that towards the origin the potential is 'heightwise truncated' (in full 2-space the potential 
is not differentiable at the origin). Thus, to harmonically approximate the potential we 
must assure that we are sufficiently far from the origin, i.e. loosely AVq must be large. 
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To obtain approximate solutions to this bent configuration equation we make the sub- 
stitution v(r) = For I ^ we may expand the effective potential about its minima 

r* = r (l + e - 3e 2 + 0(e 3 )) (A.5) 
where e = -hii} 2 — \)\ and our large r condition is refined to 

274> 2 -^(P-\). (A.6) 

Moving the left boundary condition from r = to r = — oo only introduces 0(e) corrections 
so we may obtain a ID SHO equation: 

- ^u(r) + (A + A 2 (r - r*) 2 ) u(r) = E'u{r). (A.7) 

where A = r 2 A ,2 (e-e 2 +0(e 3 )) and A 2 = A ,2 (l+3e-12e 2 +0(e 3 )). Note that the I dependence 
of the solution is implicit to the behavior of A. We find that the actual wavefunction for 
the bent configuration is: 

Ai e-^ 2 r 
»W r « ; * n i —kT H ^ V A r - r* . (A.8) 
( v /?2 n ''n r !)2 V r 

Although this solution is singular at the origin this is of little concern since the metric 
contains a factor of r (the real objects of interest are y/rv(r), of which we have a good 
approximation) . 

A. 3 Franck- Condon Factors 

We wish to calculate the overlap of wavefunctions from a linear and bent configuration. 
Given the SO (2) symmetry of both configurations the angular part of the wavefunctions 
simply gives the 5ij> selection rule. Hence, the integral we wish to evaluate is 

InrK = M <,i N n r ,i r dr r\ l \ + h-^ r2 L\l(Xr 2 )H n r(y/l(r - r*)), (A.9) 

" 

where we introduced the short hand M and for the normalizations. 

We remove all dimensionful quantities by making the change of variable u = y/Xr. 
Additionally we decide to work in the (realistic) domain where A ~ A' by setting h- = 1 + 5 
where 5 is small. Using expressions from || for the expansion 

H a (a(x-x ))= t ntS (a, - x )H n (x), (A.10) 

0<n<s 

assuming 5 and e are of the same order, and ignoring quadratic terms: 

_ 2 <L~£) r^i^^y (A.11) 
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Not surprisingly the dominant term is independent of 5 (reflecting the fact that FC integrals 
tend to be more sensitive to changes in displacement than to dilatation). Evaluating the 



integral [15| and simplifying: 




^■i iM^w . 1 { ~ roVX ' ] '" + -' r '^ + -»- (A - 12) 

B Schrodinger and Algebraic Quantum Numbers 

In the preceding sections for each I = 0, ±1,±2... we have that n r = 0,1,2,... with 
the energy relation E oc |/| + 2n r . For the u(2) chain of the u(3) model with hamiltonian 
H = Eo + kn (the algebraic model most similar to the 2D SHO) we have as a spectra E oc n. 
The branching rules imply that for each I = 0, ±1, ±2 ... we have n = |/|, |/| + 2, |/| + 4, . . .. 
Note that both the expressions for the spectra and the rules for the ranges of the quantum 
numbers agree if we make the identification: 

n — \l\ . 
n r =^±. (B.l) 

Thus the wavefunction with lowest energy within a subspace of given I is labelled by n r = 
in the Schrodinger picture and n — \l\ in the algebraic picture. 



C Scale Changes 



Our equation 2.14 differs from the results of 0] by the substitution |/| — > n. This difference 



occurs because J7J did not compute overlaps for higher energy wavefunctions since its ap- 
proximation was most valid in the low energy limit. However, their requantization technique 
works exactly for all levels of the U(2) chain when 6 = 0. We therefore may expect these 
higher energy wavefunctions to be more reliable and repeat their procedure to calculate the 
scale change for all levels in this larger domain. Proceeding to do so one finds that the 
linear \l\ dependance of equation |2.14| was an artifact of working in the ground state of each 
|/| subspace and the true correspondence has linear contributions appearing in n. 

D Schrodinger FC Factors for dilatated r* = Har- 
monic Potentials 

Using the solution to the radial part of the Schrodinger equation for the 2D SHO stated in 
Appendix |A| the FC overlaps for a dilatated SHO become: 



n r \n' r \ (AA')I'I- 



\ (\l\+n r )\(\l\+n' r )\Jo 



r r drr 2 ^e-^ r2 L^(\r 2 )Ll[(\'r 2 ). (D.l) 

J 
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The integral may be evaluated in terms of a hypergeometric function JTH] and simplified to: 



(\\')\^n r in' r \ , ynr f 2 V" +1 (D2) 



nrKA ' } \| (n r + \l\)\(n' r + \l\)\' ^ VA + A' 

a(n r . 



minK,<) (_W + n , + m _ x, / A _ A /\ n r +<-2m 



m=0 



m! (n r — m)! (n'j, — m)! yA + A' 



where we rearranged the sum to make the formula slightly more amenable to computer 
implementation. With some minor algebra and judicious manipulation of Pochhammers 
one sees that the expression has the correct A' — > A limit of 

Given the relationships of Appendix |B| and that A = a 2 the result needed for Section O 

is: 

W£) = J n SH °»-"i „_»H«i ,(«V 2 ) (D-3) 

a n r — — ,n r — — 5 — ,t 
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SHO (Fig. 2) 


SHO (Fig. 3) 


PT (Fig. 3) 


a 2 


(cm 2 ) 


1.418 x 10 1S 


1.418 x 10 18 


1.418 x 10 18 


' bent 


(cm) 


2.282 x 1CT 9 


5.459 x 10-9 


5.459 x 10~ 9 


a 2 
ff3 


(unitless) 






312. 



Table 1: Parameters for Figures 2 and 3. The parameters of Figure 2 were chosen to reflect 
the 1 B 2 — > transitions of C^. The parameters of Figure 3 correspond to a similar 
molecule with a greater bend. 
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r (arb. units) 



Figure 1: The geometry and radial potentials (of the central atom) of a triatomic model. 
In the linear configuration (a) the potential has a true minima at r = (the Hessian is 
positive definite) and there is a normal mode for each of the two degrees of freedom. In 
the bent configuration (b) the potential has a continuous manifold of minima along a circle 
surrounding the origin (the Hessian has one positive and one eigenvalue). There is only 
one non-spurious normal mode (radial vibrations or bending modes). The second mode is 
lost to the overall rotational degree of freedom. 
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Figure 2: The FC factors from a bent to linear configuration calculated numerically assuming 
the potential is harmonic about the minima. The parameters are chosen for the l B 2 — > 
transition of CS 2 (see Table 2). The factors are plotted for a fixed bent configuration ( 1 B 2 ) 
state as a function of the radial quantum number of the linear configuration for two 
subspaces of constant I. 
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Figure 3: The FC factors from a bent to linear configuration calculated for a fictitious 
molecule having similar harmonic distance scales as CS2 but a much larger bend (see Ta- 
ble 2). The factors are plotted for a potential of both harmonic (SHO) and Poschl- Teller 
(PT) type. The depth of the PT potential was chosen so that the graph shows states up to 
nearly the disassociation energy. 
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Figure 4: The FC factors from a linear to bent configuration calculated for three different 
radial displacements assuming the transition corresponds to a (1) U(3) change of basis 
within the irrep. [N] = 120; (2) U(3) change of basis within the irrep. [N] = 80; (3) 
harmonic Schrodinger potential. 
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